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In perturbative consideration of the Yang-Mills gradient flow, it is useful to introduce 
a gauge non-covariant term (“gauge-fixing term”) to the flow equation that gives rise to 
a Gaussian damping factor also for gauge degrees of freedom. In the present paper, we 
consider a modified form of the gauge-fixing term that manifestly preserves covariance 
under the background gauge transformation. It is shown that our gauge-fixing term 
does not affect gauge-invariant quantities as the conventional gauge-fixing term. The 
formulation thus allows a background gauge covariant perturbative expansion of the flow 
equation that provides, in particular, a very efficient computational method of expansion 
coefficients in the small flow time expansion. The formulation can be generalized to 
systems containing fermions. 
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1. Introduction 

As a novel method to define renormalized quantities, the Yang-Mills gradient flow 0,8 and 
its extension to the fermion field 8 have attracted much attention in recent years, mainly 

is a recent review, and Refs. [5-26] are 


in the context of lattice gauge theory. Reference 
more recent related studies. 

Although the gradient flow in lattice gauge theory is utilized to study non-perturbative 
dynamics of gauge theory, information available through perturbative theory is always useful 
because the latter is well under analytic control. In the present paper, aiming at possible 
simplification in perturbative calculations associated with the gradient flow, we consider 
the application of the idea of the background Held method [271—31] to the gradient flow. It 
is well known that this method considerably simplifies perturbative computation of, e.g., 
renormalization constants. 

As clarified in Ref. 0 ]> f° r perturbative consideration of the gradient flow, it is useful 
to introduce a “gauge-fixing term” that breaks gauge covariance of the flow equation; 
this term gives rise to a Gaussian damping factor also for gauge degrees of freedom and 
ensures a convergence property of momentum integrals. Here, we consider a modified form 
of the gauge-fixing term in the flow equation that manifestly preserves covariance under the 
background gauge transformation. It is shown that our gauge-fixing term does not affect 
gauge-invariant quantities, as the conventional gauge-fixing term. This formulation thus 
allows a background gauge covariant perturbative expansion of the flow equation that pro¬ 
vides, in particular, a very efficient computational method of expansion coefficients in the 
small flow time expansion [ 2 ]. 

This paper is organized as follows. In Sect. El we present our general formulation. Both 
flow equations for the gauge field and for the fermion fields are considered. The most impor¬ 
tant observation is the independence of gauge-invariant quantities on the gauge-fixing term 
we introduce (Sect. 12.31) . In subsequent sections, we consider applications of the formula¬ 
tion: In Sect. E we consider the computation of expansion coefficients in the small flow 
time expansion 2] relevant to the construction of the lattice energy-momentum tensor; this 
computation was carried out in Ref. 32fl using a cumbersome diagrammatic method. We 
observe that the application of our formulation provides a very efficient non-diagrammatic 
computational method, that is quite analogous to that of Ref. [33|, for the expansion coef¬ 
ficients Q In Sect. El we consider the small flow time expansion relevant to the construction 
of the axial-vector current 34j. The last section is devoted to the conclusion. 

Here is a summary of our notation: Our generators T a of the gauge group G are anti- 
Hermitian and the structure constants are defined by [T a ,T b ] = f abc T c . Quadratic Casimirs 
are defined by f acd f bcd = C2(G)5 ab and, for a representation R , trf i(T a T b ) = —T(R)5 ab 
and T a T a = — We also denote tr#(l) = dim(R). For example, for the fundamental 

N representation of SU(N ) for which dim(IV) = N, the conventional choice is 


C 2 {SU{N )) = N, T(N) = 


1 


c 2 (N) = 


N 1 2 - 1 
2 N ' 


( 1 . 1 ) 


1 Unfortunately, the, results of this new simple computational scheme do not coincide with 

the results in Ref. 32], revealing that there are errors in the one-loop diagrammatic calculation 
in Ref. [321 ]. The diagrams in which the mistakes wer e made in Ref. [32] have been completely 
identified. For corrected results, see the errata for Refs. [32|, [35| and Refs. 136l I37|. 
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Our gamma matrices are Hermitian and for the trace over the spinor index we set tr(l) = 4 
for any spacetime dimension D. The chiral matrix is defined by 75 = 70717273 for any D 
and thus 


t r (757M7i'7p7er) 


4 e^pa, p,v,p,a <E {0,1,2, 3}, 
0, otherwise, 


( 1 . 2 ) 


where the totally anti-symmetric tensor is normalized as €0123 = 1- 


2. Flow equations with a background covariant gauge 

2.1. Gradient flow equation with a background covariant gauge 
The gradient flow for the gauge potential is defined by [lj 

dtBflt, x ) = D u G v flt, x ) + a 0 D^d u B u (t, x), Bflt = 0, x) = Aflx), (2.1) 

where t € [0,oo), and 

Du = df, + [Bft,-], G^(t,x) = d^Bflflx) - d u Bflt,x) + [Bflt,x),B u (t,x)] (2.2) 


denote the covariant derivative and the field strength of the flowed gauge field, respectively. 
The last term in the first relation of Eq. (12.11) breaks gauge covariance and here it is referred 
to as a “gauge-fixing term”. As noted in Ref [if] , for perturbative consideration of the gradient 
flow, such as that in Ref. 0], it is useful to introduce such a gauge-breaking term because 
it gives rise to a Gaussian damping factor also for gauge degrees of freedom and ensures a 
convergence property of momentum integrals. It can, however, be shown that [lj any gauge- 
invariant quantity, that does not contain the flow time derivative dt , is independent of the 
“gauge parameter” ao and physical observables are not affected by the gauge-fixing term. 

In the present paper, we propose a slight modification of the gauge-fixing term in Eq. (12.11) . 
First, following the general idea of the background held method [27431]], we decompose the 
original gauge potential into the background part Aflx) and the quantum part aflx) as 

Aflx) = Aflx) + aflx). (2.3) 

We also decompose the flowed gauge potential B(t,x ) into the background part Bflt,x) and 
the quantum part b^(t,x) as 

Bflt,x) = Bflt,x)+ bflt,x). (2.4) 

Then, our proposal is to adopt, instead of Eq. m, 

d t Bflt, x) = D v G v flt, x) + aoD^D v b v (t, x), Bflt = 0, x) = Aflx), (2.5) 


where 

Dn = + [B^, ■] ( 2 . 6 ) 

are the covariant derivatives with respect to the background flowed field. 

As a further natural assumption, we suppose that the background held is evolved by its 
own how equation: 

d t Bflt, x) = D u G v flt, x), Bflt = 0, x) = Aflx), (2.7) 


where 

G^{t,x) = dfj,B v (t,x) - d u Bflt,x) + [Bflt,x),B u (t,x)\ (2.8) 

is the held strength of the background held. 
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2.2. Covariance under the background gauge transformation 
The original gauge transformation 


A^x) ->• Afj,(x) + DnU}(x), 


(2.9) 


may be decomposed into the background part and the quantum part; how this decomposition 
is made is the heart of the background field method |27H3l|. A fundamental notion is the 
background gauge transformation , defined by 


A^{x) ->• Afj,(x) + Dfj,uj(x), a^(x) ->• a M (a?) + [a fl (x),uj(x)}. 


( 2 . 10 ) 


The sum of these two reproduces the original gauge transformation (ZSD- Under this 
background gauge transformation, the quantum gauge field transforms as the adjoint 
representation. This transformation can also be generalized to the flowed fields aJ! 


Bn(t,x) ->• BfiftiX) + b,j,u(x), bfj,(t,x) bfj(t,x) + [b fl (t,x),uj(x)\. (2.11) 


Note that here we are assuming that the transformation function uj(x) does not depend on 
the flow time t. 

Since D M in Eq. (12.61) transforms covariantly under the background gauge transforma¬ 
tion (12.111) . our flow equation (12.51) transforms covariantly under the background gauge 
transformation; fields transformed by the background gauge transformation obey the identi¬ 
cal equation. This is the reason for our choice of the particular gauge-fixing term in Eq. (|2.5I) 
instead of the conventional one in Eq. (12.11) . 

Now let us confirm that our gauge-fixing term in Eq. (12.51) does not affect gauge-invariant 
quantities. 


2.3. Independence of gauge-invariant quantities of the gauge parameter ao 
Although our gauge-fixing term aoD^D u b u (j:,x) in Eq. (12.51) differs from the conventional 
one in Eq. (12.11) . one can still see that any gauge-invariant quantity, that does not contain 
the flow time derivative dt is independent of the “gauge parameter” ao; the gauge-fixing 
term thus does not affect gauge-invariant quantities. 

To see this, we consider the following “quantum gauge transformation”^ 

Bfj,(t,x) —>• Bfj(t,x), bfj,(t, x) —>■ bp(t, x) + D^ui(t, x), (2.12) 

whose transformation function uj(t,x) does depend on the flow time t. Note that the sum 
of these two reproduces the original gauge transformation (I2.9jl with uj(x) —>• u(t,x). Under 


2 The covariant derivative D^ in the first relation is defined with respect to the flowed background 
field x). 

3 Here the covariant derivative D^ in the second expression is defined with respect to the flowed 
field B^(t, x ). 
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this infinitesimal transformation, we have 

d t B^(t, x) d t Bflt,x) + [d t B^t,x),uj(t,x)} + D^dtufl, x), (2.13) 

D u G Ufi (t,x) ->■ D u G Ufl (t,x) + [D u G UIJ ,(t,x),uj(t,x)], (2.14) 

D u b u (t, x) ->• D p b u (t, x) + D u D u uj(t , x), (2.15) 

and 

D^D v b v (t, x) ->• D^ L D v b u {t 1 x) + [D^D u b u (t, x),u(t, x)] + D^D u D u u(t, x), (2.16) 

where in deriving the last expression we have noted 

D u D„uj(t,x) = D u D u uj(t,x) + [D u b u (t,x),oj(t,x)\. (2.17) 

From these expressions, we see that under Eq. 112.121) . the flow equation (| 2.5 1) changes to 

dtBflt, x) = D u G vll {t, x) + aoD^D v b v (t, x) - D^(d t - a 0 D u D u )uj(t, x). (2.18) 

This shows that, by choosing u(t,x) as a solution of 

{dt — aoD u D u )u(t, x) = — 5aoD u b u (t, x), u(t = 0,x) = 0, (2-19) 


the transformed flowed field (that has the same initial value as the original one) obeys the 
flow equation (12.51) with «o —» ao + £ao- Since a gauge-invariant quantity that does not 
contain the t derivative is invariant under Eq. (12.121) , this shows that such a gauge-invariant 
quantity is independent of ao- Physical observables are not affected by the gauge fixing 
term in Eq. (12.51) ; the introduction of the gauge-fixing term is thus a physically allowed 
modification of the flow equation. 


2-4- Classical perturbative solution to the flow equation 

Now, using Eq. (12.71) in Eq. (12.51) . we have the flow equation for the quantum field: 


d t bflt,x)= d^D 2 + (a 0 - l)D^D v b v (t, x) + 2[G /iI/ (t, x), b v (t, x)] + R^t, x), (2.20) 

where 


Rflt,x) = 2 [b v (t,x),D v bn(t,x)\ - [b u (t,x),b^b u (t,x)\ 

+ (a 0 - l)[bflt,x),D u b v (t,x)\ + [b v (t,x), [b„(t,x),b,j,(t,x)]\. (2.21) 

The adjoint actions in these expressions can conveniently be expressed in terms of matrix 
multiplication, if one introduces the adjoint representation. We thus define 

Bf(t,x) = B^t,x)f ac \ (2.22) 

Vf^5 ab d ll + Bf{t,x), (2.23) 

gfi(t,x) = G^(t,x)r b . (2.24) 


With these notations, the flow equation for the quantum field reads 


I ab 


d t b a (t,x)= S^V 2 + (a 0 - l)VflD u + 2Q flv (t,x) b b u (t, x) + R a (t, x), 


(2.25) 
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where 


R a Jt,x) = 2 f abc bl(t,x)Vt d b d M,x) - f abc b b u (t,x)V c X(t,x) 


cabcib , 


+ (ao - 1 )f Mc b b li (t,x)V^bi(t,x) + f^n^xK^x^x). 


(2.26) 


A formal solution to Eq. (12.2511 is then given by 


K(t,x) = / d 4 y 




<26 


(2.27) 


where the heat kernel K" b (x, y)^ u is defined as an object that satisfies 


d t Kf(x, y)^ = S /lX V z + (a 0 - \)VfiD x + 2 x) K?(x, y) Xu 


>cb( 


K£ 0 (x,y)^ = 6 ab 6 lx „8(x-y). 


(2.28) 

(2.29) 


The heat kernel defined by Eqs. (I2.28P and (12.291) may be expressed in the form of a time- 
ordered product containing the flowed background field. If one is considering a particular 
situation in which the background field A(x) can be assumed to obey the Yang-Mills equation 
of motion, 

D u F Ufl (x) = 0, (2.30) 

then Eq. ()2.7j) implies that the background gauge field does not flow: 

B(t,x) = A(x). (2.31) 


Then the heat kernel in the “Feynman gauge” ao = 1 can be written, suppressing the gauge 
and Lorentz indices, in a very compact form: 

r rt 'i 


K t (x,y) = T exp 


ds 


i>l + 2 G(s, x) 


S(x - y) 


= e t[Dl+2f( x )\ 6{x _ y) 


(2.32) 


In the last expression, the covariant derivative is defined with respect to the background 
gauge field at vanishing flow time, A /t (x); we have also introduced the corresponding field 
strength in the adjoint representation, 

^(x) = (2.33) 

In the application to the small flow time expansion in the next section, we can assume 
Eq. (12.301) without loss of generality. We can then use Eq. (12.32j) for the heat kernel which 
greatly simplifies the computation. 


2.5. Tree-level propagator of the flowed gauge field 

So far we have considered the flow equation (12.201) at the classical level. The quantum field 
at vanishing flow time, a b (y ), contained in Eq. (12.271) is actually the subject of the functional 
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integral with the Boltzmann weight, specified by the Yang-Mills action 

1 


S = 




d (x)F“ (x) 


and the gauge-fixing term in the background gauge 27H31I] 


5, 


gauge-fixing 


Ao. f 

H J ' 


d 4 x D^a a flx)D u a^(x), 


(2.34) 


(2.35) 


which also preserves covariance under the background gauge transformation. Then, in the 
presence of the background field, the action quadratic in the quantum field is given by 


(S + ^gauge-fixing ) 1 0 (a 2 




ab 


<Txa a Jx) Sfjv'D 2 + (A 0 - l)V^D u + 2F llv {x) <(x), 


and thus the tree-level propagator in the Feynman gauge Aq = 1 is written as 


<(zK(y )) 0 = so 


-l 


Vl + 2F{x) 


5(x - y ). 


(2.36) 


(2.37) 


pis 


If one can further assume Eq. f)2.30ji for the background field, then the heat kernel (in the 
Feynman gauge) is given by Eq. (12.321) . Then, from Eq. (12.271) . the tree-level propagator of 
the flowed quantum field, in the presence of the background field, is given by 


b;(t,x)b b u (s,y)) a = gl [e^ 2 ^ 


-1 


Vl + 2F{x 

{t+s)[T>l+ 2 T(x)] 


/"e s P^+ 2 -^hd]'j d (x — y) 

\ J ISp 

pp 


= 9 o 

where in the last equality we have noted 

'-> 2 y + 2F(y)] ab 5{x- 

J pis L J isp 

The above expression (I2.38p . which is manifestly covariant under the background gauge 
transformation, will be fully employed in our application to the small flow time expansion 
in the next section. 


-1 A 

ab 

5{x-y), 

pis 

(2.38) 

Vl + 2F{x)J 

V 2 x + 2 F(x) 

ba 

5{x - y). 

(2.39) 


2.6. Fermion flow 

We can also consider the “background covariant gauge” in the flow of fermion fields Qj; we 
adopt the following flow equations: 

dtx(t,x) = {z? 2 - a 0 [£>PPt,x)]} x(t = 0, x) = V>(x), (2.40) 

dtx(t,x) = x(t,x) j^> 2 +a 0 [Dppt,x)]} , x(t = 0, x) = $(x), (2.41) 

where the covariant derivatives on the fermion fields are defined by 

D fl = d /t + B t j, ^ — (2.42) 

and 

Dp = < 9 ^ + B^, ^5^= d^ — B^. ( 2 . 43 ) 

On the other hand, in these expressions and in what follows, [Dp>^(t, x)\ stands for the 
background covariant derivative on the quantum gauge field, defined in Eq. (12.61) . 
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One can again see that the gauge parameter «o is irrelevant for gauge-invariant quantities 
with our gauge-fixing terms in Eqs. (12.401) and (12.411) . To see this, we again consider the 
infinitesimal transformation, Eq. (12.12|1 . and 

X(t,x) -> [1 -u(t,x)\x(t,x), x(t,x) -> x(t,x) [1 + u(t,x)] . (2.44) 

Then using Eq. (12.151) . after some calculation, we see that the flow equations are changed as 
dtx(t,x) = [d 2 - a 0 [D^b^(t,x)]j x(t,x) + (d t - a 0 D /J ,D ll )uj(t,x)x{t,x), (2.45) 

dtx(t,x) = x(t,x) j^ 2 + oco[D tl b IJ ,(t,x)]j - x(t, x)(d t - aoD^D^u}^, x). (2.46) 

These show that again by choosing u(t,x) as the solution to Eq. (12.191) . we can shift ao 
to «o + boio- Gauge-invariant quantities (that do not contain the flow time derivative) are 
hence not affected by the gauge-fixing terms in Eqs. (12.40J) and (12.41 J) . 

We also decompose the fermion fields into the background part and the quantum part as 


X(t,x) = x(t,x) + k(t,x), 

X{t,x) = x(t,x) + k{t,x). 

(2.47) 

lf>(x) = lj)(x) +p(x), 

'ijj(x) = ijj(x) + p{x), 

(2.48) 

and assume that the background fields themselves 

are evolved according to 


dtx(t,x) = D 2 x(t,x), 

><> 

c-+* 

II 

o 

II 

(5) 

(2.49) 

dtX{t> x ) = X(t,x)i5 , 

X(t = 0,x) = ijj(x). 

(2.50) 


Then, from Eqs. (|2.40|) and (|2.4ip . the quantum fields obey the flow equations 
d t k(t,x ) = | D 2 - a 0 [l) M & M (t,:r)]} k(t,x) 

+ {(1 - ao^Dpbftfax)] + 2b tl (t,x)D fl + b^t, x)b^(t, x)| x(t,x), 


k(t = 0, x) = p(x), (2.51) 

d t k(t,x) = k(t,x ) ^D 2 + a 0 [D^bfj,(t,x )]} 

+ X(t,x) I -(1 - a 0 )[T) M 6 /i (i,a:)] - 2D ft b tl (t,x) + & M (i, x)^^,®) i , 


k(t = 0, x) = p(x). (2.52) 

If we further assume that the background gauge field fulfills Eq. (12.301) . the background 
gauge field does not evolve as Eq. (12.311) and we can write down relatively simple expressions 
for the solution of the fermion flow. The solution to the flow equations (12.491) and (I2.5UD can 
be expressed as 

X(t,x) = e tD2 ^{x), x{t,x) = i>{x)(fi . (2.53) 

Then, the solution to the flow equations (12.511) and (12.521) is given by 


k(t, x) = e t22>2 p(x) + / dse^ s ^ 2 2b fi (s, x)D^ + b^s, x)b^(s, x) e s ^ 2 %J>(x) + k(s,x) 


k(t, x) = p(x)e t< ^ + [ ds 


il)(x)e sD +k(s,x) 


2D^bJs, x) + 6 m (s, x)6 m (s, x) 


(2.54) 

(2.55) 
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where we have adopted the “Feynman gauge” op = 1 for simplicity. 

The quantum fields at vanishing flow time, p(x) and p(x), are subjects of the functional 
integral with the conventional action, 


S = J d D x'il>(x)(]p + mo)i/>(x) 

= j d D x ip(x)+p(x) (Jj) + c/l-\- mo'j $(x) + p(x) 


(2.56) 


Thus the tree-level propagator of quantum fermion fields, in the presence of the background 
gauge field, is given by 

1 


(p(x)p(y)) o = 


px + rn 0 


S(x - y ). 


(2.57) 


3. Application: Small flow time expansion relevant to the energy-momentum 
tensor 

As noted in Ref. [2], any local composite operator of flowed fields can be expressed as, in 
the limit of t —>• 0, an asymptotic series of local composite operators of fields at vanishing 
flow time. In Ref. [33], this small flow time expansion (with use of perturbation theory) was 


exploited to construct a universal formula for the energy-momentum tensor. This formula 
with lattice regularization was then numerically tested in Ref. [38] by applying it to the bulk 
thermodynamics of quenched QCD. The universal formula can be generalized to general 
vector-like gauge theories [35] and to various asymptotically free theories [39 4lj]. In Refs. 


42 


431 ], application of the gradient flow to the lattice energy-momentum tensor is studied from 
a somewhat different perspective. 


In the present paper, we restrict ourselves to the case of the pure Yang-Mills theory 32] 
and consider the small flow time expansion in the forinQ 


G a iip {t 1 x)G a up (t,x) 

<G^ p (t,a:)G^ p (t,®)) + Cn(t)^ p (a:)^p(a:) + + O(t), (3.1) 

where the 0(t ) term is the contribution of composite operators of the mass dimension being 
equal to or greater than six. Because of symmetry, only the above two four-dimensional 
operators can appear on the right-hand side. The expansion coefficients can be evaluated in 
perturbation theory and, to the one-loop order, we write 

Cn(t) = 1 + (n(t) + • • • , C 12 (t) = 0 + (^(t) + • • • • (3-2) 

If these coefficients are obtained in the dimensional regularization (with the spacetime dimen¬ 
sion D = 4 — 2e), then the correctly normalized conserved energy-momentum tensor (with 


4 In Appendix [B] we compute the small flow time expansion of an operator corresponding to the 
topological density—another gauge-invariant dimension-four operator. 
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the vacuum expectation value subtracted) can be written as 32|, l35j] 


{T, w \r{x) = lirrij d (t) 


t—>o l 


pya 


(t,x)G“ p (t,x) - -SpuG pa (t, x)G a pa (t,x) 


+ c 2 (t) [S^G a pa (t,x)G a pa (t,x) - (5 pu G a pa (t,x)G a pa (t,x))] , (3.3) 


where 


c i(0 = 4 l 1 ~ Cn } W 

So L 


1 


-{k'So) 


(3.4) 


c 2 (t) = -2 
So 

Since bare composite operators of the flowed gauge field are automatically renormalized 
operators Q], the formula (13.31) should hold irrespective of an adopted regularization J§ in 
this sense the formula is universal. In particular, it should hold with lattice regularization 
with which the construction of a correctly normalized conserved energy-momentum tensor 
is not straightforward. It is thus of great interest to compute the expansion coefficients 
in Eq. (13.11) . As we will see below, the background field method we have developed provides 
a very efficient non-diagrammatic computational method of the expansion coefficients (at 
least in the one-loop level). 

Now, to determine the expansion coefficients £n(t) and in Eq- (j3.ll) . we consider 

1PI diagrams containing the composite operators G pp (t, x)G^ p (t, x) and F pp (x)F£ p (x) with 
external lines of the background gauge field B p {t , a;) only (i.e., no external line of the quantum 
field). In the tree level, since the flow time evolution is purely classical, 


<G“ p (t,s)G“ p (t,x)) lpi f; p (*)F“ p (z) + 0(t), 
<e; p (x)e“ p (x)) 1p1 = f; p {x)fz p {x). 


(3.5) 

(3.6) 


Comparing these two relations, we find 

Gp P (t,x)G* p (t,x) *2° F^(x)F^(x) + Q(t). (3.7) 


This gives the tree-level contributions in Eq. (13.21) . 

Next, we consider one-loop 1PI diagrams containing the composite operators with exter¬ 
nal lines of the background gauge field. Such 1PI diagrams can be obtained, by taking 
the contraction of quantum fields in the expansion of the composite operators by the 
propagator in the presence of the background field. The expansion of the composite 
operator G pp (t,x)G® p (t,x) in the quadratic order yields 


G £p( i >z)G“ p (t,x)| 0(b2) 

r , -| a r A 

— {^fia^uS^/3'y Yr) G\~/ T ^F a bfj(t,x') F ft A,,(t, X) 

- b u (t,x)J r p , p (x)b p (t 1 x) - b p {t : x)F up {x)b p {t,x). 


(3.8) 


At this stage we note that to read off the expansion coefficients in Eq. (13.11) as Eq. (13.71) . we 
can assume that the background field satisfies the Yang-Mills equation of motion, Eq (12.301) , 
because F pp (x)F£ p (x) does not vanish under the equation of motion. Then the background 


5 The coefficients c\(t) and c 2 (f) in Eq. (13.411 become finite for 
renormalized quantities; see below. 


0 when expressed in terms of 


10 



















field does not flow B(t,x ) = A(x) and we can use the simple expression (|2.38l) for the 
propagator. The contraction then yields 


x )G up (t, x)|o(fi 2 )/ 1PI 
— <7o( ( W ( W ( Vy T S 


ab I n 2tA ^ \ 


be 


T V?5{x-y) I 


,2tA 

a; 


ba 


+ g^ p [x) ( e^- ) 6(x - y)\ v=x + g^ p (x) ( e 


pv 


2tA 


A 


A J 

ba 


Pi 


y=x 


5{x-y)V 


y=x i 


(3.9) 


pp, 


where 


A = V 2 + 2 F, 


(3.10) 


and we have noted 


V%8(x-y) = -V%8{x-y). 


(3.11) 


However, as explained in Refs. [37] (see also Ref. [33]), instead of the 1PI function (13.91) 
itself, it is much convenient to consider the difference 


(G pp (t,x)G“ p (t,x)\ 0{b2) 



(3.12) 


because possible infrared divergences are cancelled out in this combination. The one-loop 
1PI function, which contains F pp (x)F£ p (x) is given by simply setting t = 0 in Eq. (13.91) . Then 
the difference can be expressed as an integral over an auxiliary variable £ as, 


(G pp (t,x)G“ p (t,x) | 0(ft2) 


= 2<?o 





Fp P {x)F vp (x) | 0(a 2 )) ipi 

fipadu-ySps — dppdvsda'Y + ^p/3^u'yS a s)V^ > (e 2 ^') Vf l 

\ J p'y 


+ ^(x) 


3 2 ? A\ 


ba 


pv 


+ 




ba 


pp 


S(x-y) L =; 


y=x 


(3.13) 


In this expression, it is obvious that there is no infrared divergence, because derivative 
operators appear only in positive powers. 

We then set 


S(x-y) = f (3.14) 

and moves the plain wave e ipx the most left-hand side, by noting 

V p e ipx =e ipx (i Pp + V p ), 


11 


(3.15) 









as usual in the calculation of anomalies in the path integral [33), I44M471]. After the rescaling 
of the integration variables, —>• Pp,/y/£, we have 


G a pp {t,x)G“ p (t,x)\ Q(b2) - FZ p (x)FZ p {x )| g(o2)/ ipi 


29? /‘d«- D/2 / 


d D p 


/ o 

x tr 


(2vr) 




2p 2 


&(ia&i' r y&/3S $p.fi^uS^a'y 4“ 3p,/3^i''y&oi5 ) 

x C 1 ( iPa + V^a) ^Vtp-v+2t i A.\ ^ + y/lv^j 


+ F^ p (x) [e 


3 4iV?p-i>+2?A4 


j +J^ p (x)(^e' 


4i7Cp-X)+25A^ 


Pfi 


(3.16) 


where the trace tr is for the gauge indices. It is interesting to note that all the information 
of relevant one-loop IPI diagrams is contained in this single compact expression; in a con¬ 
ventional calculational scheme [32j], on the other hand, one has to compute at least 12 IPI 


diagrams to obtain the expansion coefficients in Eq. m- 
Now, since we are interested in the small flow time limit t —» 0 of Eq. (13.1611 and since £ € 
[0,f], we may expand the integrand with respect to £. For t —» 0, only terms to 0(^~ D ^ 2+1 ) 
under the integral can give rise to non-vanishing contributions for D —> 4. The expansion 
of the combination e 4 * v ^ p '- E>+2 ^ A to 0(£ 2 ) is given in Appendix [XJ Although the remaining 
algebraic calculation after the Gaussian integration over p, by noting 


[D p ,T> a \ = F f 


per 


(3.17) 


is somewhat lengthy, it is rather straightforward. In this calculation, it is quite helpful to 
note that the final expression for Eq. (13.1611 must be symmetric under p o v by definition; 
we may thus simply discard any terms anti-symmetric under p <->• v. In this way, we finally 
arrive at 


Gp P (t,x)G“ p {t,x)\ 0{b2) - F“ p (x)F? p (x)\ Q(a 


0(a 2 ) 


IPI 




+ 


So 


(4tt) 


+ 


11 , , i 7 
— e(t _1 + - 
3 W 3 


tr 


Hxf 


J /IV 


So 


(4?r) 


11 / \_i 1 

- e(t) - 

12 w 6 


Spu tr 


F{x) 2 +0(t), 

J pp 


(3.18) 


where 

Since 

tr 


= c 2 (g)f“ (*)F“ (x), (3.20) 


e(t) * =— |-ln(87rt). 


each rbda r?c 


(3.19) 


J flU 
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recalling the tree-level relations ([3.51) and (13.6p . Eq. (|3.18|) shows that the small flow time 
expansion to the one-loop order is given by 


G^ P {t,x)Gt P {t,x) 


up\ 


t-> o 


^dim (G)^V 


+ 1 + 


9o 


(47r)' 


;C 2 (G) 


+ tA^C 2 (G) 


(Ait)'* 


11 , . i 7 
3 W 3 

11 , \_i 1 

- 12 eW ' “ 6 


f; p (x) f 2 p (x) 
5, u F a (x)F a (x) + 0(t). 


From this, £ j| and hr Ed- (13.21) are given by 


di) : 


dR) = (Sy C2(G) 


dl(t) = 


9o 

(4vr) 


■C 2 (G) 


11 , x-i 7' 

Y e(t) + 3 

11 / 3-1 1 

-l2 e(t> - 6 


(3.21) 

(3.22) 

(3.23) 


and then Eq. (13.41) gives ci(t) and c 2 (t) in Eq. (13.31) . In terms of the renormalized gauge 
coupling g in the MS scheme, 


1 


1 


1 


— = -y + &o-In /i 

9o 9 \ e 

where bo is the one-loop coefficient in the beta function, 

6 » = (Yy C 2 (G) f' 


we have 


ci(t) = -j~ b 0 ln(8vr/u 2 t) - jj-y 2 C 2 (G) 1 -, 
r (4=71 -y 3 

c 2 (£) = ^0- 


(3.24) 

(3.25) 

(3.26) 

(3.27) 


The above non-diagrammatic one-loop computation of coefficients ci(£) and c 2 (t) is much 
simpler and quicker than the diagrammatic calculation carried out in Ref. [32j. Unfortu¬ 
nately, the results of the above calculation do not coincide with the results in Ref. [321 ]. 
revealing that there are errors in the one-loop diagrammatic calculation in Ref. 32] § 


6 In particular, Eq. (4.30) and (4.31) of Ref. ]32| should be 


Cl = lny/n+ — ~ 0.890547, 

22 

c 2 = ln ^-^7 + 4^ 0-834762. 


(3.28) 

(3.29) 


Equations (4.32) and (4.33) of Ref. 35 should be replaced by Eqs. (13.221) and (13.231) . respectively, 
and consequently, Eq. (4.72) of Ref. [35| should be 


ci (t) 


1 

g(i/V8tr 


bo In 7r — 


1 

(4tt)2 


\c2(G) 


\T(R)Nt 


(3.30) 
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4. Application: Small flow time expansion of the axial-vector current 


As another application of the present formulation, we consider the small flow time expansion 
of the axial-vector current of the flowed fermion fields [341 ]: 


X(t, z) 7 m 75 t A x(t, x ) 


H-£ (1) (*)l ^(z) 7 //y 6 t A ip(x) + 0(t), 


(4.1) 


where t A is the generator of the flavor symmetry group and (t) is the expansion coefficient 
at the one-loop level. Because of symmetry, only the axial-vector current at vanishing flow 
time can appear as the leading 0(t° ) term in the right-hand side. To find £^(t), we set the 
background gauge field to zero and consider one-loop 1PI diagrams containing the composite 
operator x(f, 2 ) 7 ^ 75 t A x(t, x ) with external lines of the background fermion fields, x) 
and ip(t,x ) (no external line of the quantum fields). 

As Eqs. (|3.5I) and (|3.6lh at the tree level, 

(x(t, z)7m75* A xO, z)) 1PI + 0(t), (4.2) 

($(®)7M75t J V(®)) 1PI = ^{x)^^t A ^(x). (4.3) 

For the one-loop level, again as Eq. (13.12f) . it is convenient to consider the difference: 


(x(t, x) 7 m 75 t A x(t, x) - ^{x)^^ b t A ^{x)) im , (4.4) 

because this combination is free from infrared divergences. For the first term of Eq. (S3D, 
by using the decomposition (12.4711 . we have 


x(t-,x)7^ b t A x{t,x) 

= X(t, 3 ?) 7 m 75 t A x(t, x) + x(t, 3 j) 7 m 75 t A k(t, x) + k(t, x)x^t A x(t, x) + k(t, x)x^t A 


k(t, x). 
(4.5) 


We can then use Eqs. (12.54|) and (I2.55P to express the quantum flowed Helds, k(t, x) 
and k(t,x), in terms of fermion fields at vanishing flow time. Equation (14.11) shows that 
to find the coefficient ^\t), we may set the background fermion fields to be constant which 
makes the calculation quite easy. Then, as terms which contribute to one-loop 1PI diagrams 
with external lines of background fermion fields, we have 

k(t, x) = e td2 p(x) + f dse^" s ^“ 6 M (s,x) 6 M (s,x )-0 + f dse^^ 92 2b li (s,x)d IJ ,e s92 p(x), 

Jo Jo 

(4.6) 


k(t, x) 

= p(x)e td2 + [ ds# M (s,x) 6 M (s,a;)e {< ~ s)9 + [ ds p(x)e s 92 (-2)^ ^b^s, x)e (f ~ s) 9 

Jo Jo 


(4.7) 


where g(l/\/&i) is the running gauge coupling in the MS scheme at the renormalization scale fi = 
l/\/8 1. The expressions just below Eqs. (5) and ( 6 ) of Ref. [H| should be 

sr = ^ + \p E - In 2 ~ -0.0863575299274, (3.31) 

21 b 27 

S 2 = fr-x 79 = 77 rr - 0.0557851239669. (3.32) 

44 2b\ 484 v ’ 

The erratum for Ref. i38:j will appear soon. 
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In Eqs. (14.61) and ([4.7)1 . the quantum fields at vanishing flow time, p(x) andp(x), are subject 
to the functional integral with the action (|2.56p . Through the interaction terms in Eq. (12.561) . 
p(x) and p(x) become the background fields, il> and tjj. Considering the contraction by 
the propagator (12.571) in (p(x)(—l)f d D yp{y)(/L{y)'tp) and ((—1) f d D y '$4(y)p(y)P( x ))^ this 
effect of interaction vertices can effectively be represented by the substitutions, 

p(x) -> - [ d D y — - 5{x - y)<ji(y )^, (4.8) 

J p x + m 0 

P(x) -> - f d D y^(y )-—J- 5(y — x). (4.9) 

J <p y + nriQ 

Note that these substitutions accompany a gauge interaction vertex coming from the 
action (12.561) because we are considering 1PI diagrams. Then the contraction of the quan¬ 
tum gauge fields in the expectation value of Eq. (14.51) by the propagator (12.381) is very 
simple. In this way, we have the one-loop expression for the first term of Eq. (14.41) . Then, 
by simply setting t = 0 in that expression, we have the one-loop expression for the second 
term of Eq. (14.41) . The resulting difference is free from infrared divergences and under the 
dimensional regularization, it is simple to obtain at the one-loop level 




ipi 


t-y o g o 


(47t)' 


:C 2 (R)(-3) 


<*r L + 


'0(^)7/^ 75 t A i>(x) + 0{gl). 


Because of the tree-level relations (14.21) and (14.31) . Eq. (14.101) shows that 


t-y o 

r^j 


x(f,x) 7 At 7 5 t x{t,x) 
which coincides with the result in Ref. 


1 + (£y C2(R)( “ 3) 


e(t) + 6 


^(x)lLLlbt A ^{x) 


(4.10) 


(4.11) 


As discussed in Ref. 


34J]- 


34|, Eq. (14.111) shows that the correctly normalized axial-vector 
current can be expressed as 


4(x) = lim(l + 9{l /^ Y C 2 {R) 


•?5 fj, 


t-y 0 ( 


(4vr) 


-i+ln(432) 


X.(t,x)'Yu'y5t A x( t > x ), 


(4.12) 


where 


m 


X(t, x) = 


X(t, x ) = 

(Nf is the number of flavors) and 


-2dim(fljty- — x (t,x), 


X(t,x), 


\ (4vr) 2 t 2 (x(t, x) p x(t, x) 


—2 dim (R)Nf 


T7 - 


\ (4vr ) 2 t 2 {x{t,x)If)x(t,x ) 


*3/1 = D tl - 


(4.13) 


(4.14) 


(4.15) 


7 Note that we are now setting the background gauge field to zero. 
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5. Conclusion 


In the present paper, we have developed a background field method (or a background gauge 
covariant gauge fixing, more appropriately) for the gradient flow equations. This formu¬ 
lation allows a manifestly background gauge covariant perturbative expansion of the flow 
equations. We illustrated the power of the method by applying it to the one-loop calculation 
of expansion coefficients in the small flow time expansion relevant to the energy-momentum 
tensor. This new simple computational scheme revealed that there were errors in the old 
diagrammatic calculation in Ref. |32|] (the errors have been identified and corrected |36|,|37j). 

Since our method provides a greatly simplified computational scheme for known one-loop 
computations, we can expect that it can also be useful in more complicated situations, such 
as the two-loop computation of the expansion coefficients. We hope to come back to possible 
further applications of the present formulation in the near future. 
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A. Expansion of 

A straightforward expansion yields 

/' e 4iV&-X>+2£A\ 

V / fj,i/ 


= Sfiu + 4 ip ■ V5^ 1/2 + 2A M1/ - 8 (p ■ V) 2 5 


£ 


+ 


+ 


32 


4 iA^p -V + Aip- T>A pv - y i(p ■ V) z 5^ 


£3/2 


2A^ - a: A ^(p • Vf - a: p . vA pvP VfA^ + ^(p ■ V)*5\ 


16 


16 


16. 


32 


— pus \xr- — j tr — —\xisr ~ ^ \jr — j — f-iv > ^ \jr — j J flU 


+ o(e 5/2 ), 


r 

(Al) 


where 


A/,1, — 't) 2 bp,v + 2 J- f (A2) 
A ;= V 2 V 2 5^ + 2V 2 T pv + 2 f pv V 2 + 4 f 2 „. (A3) 


B. Small flow time expansion of the topological density 

In this appendix, we present the small flow time expansion of an operator corresponding 
to the topological charge density in the one-loop order. First, to the quadratic order in the 
quantum field, we have 


tfivpcr (t, x)G pa (t , x) \q^ 2 


— “2^pi/pa > 2 


12 V p b v (t,x) V p b a {t,x) -b l _ l (t,x)f up (x)b C7 (t,x) j. (Bl) 
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The contraction by the propagator (12.381) then yields 


~pupaG pu (t, x )G po -(t, x)\ q^, 

1 


= - 2 5 0 e 


pupa 


2f) ab I ^2£A 


A 


)/ ipi 

be 

V™ + 7^(x)[e 


up 


2tA 1 1 


A J 


ba 


per 


S(x-y)\ y=x . (B2) 


The same procedure as led to Eq. (13.161) in the main text then gives rise to 


pupaG p U (t, x)G pa (t, x)\ 0 e pupcrF fll/ (x)F pcr (x)\ 0 ^ a2 ^ j ^ , 
, rt der D/2 / d ° p c~ 2p2 

pupa I ? j (2^0 


= - 4 ffo e 


x tr 


2r x (ip + v^j 


^i^p-T>+2^A\ 


\W + 

J up 


V&>) +Fpu(x) ^e 4 ^^+ 2 « A ) 


pa 

(B3) 


The expansion with respect to £ is much simpler than Eq. (13.161) . Thus we give some details 
of the calculation for illustration. 

First, in the integrand of Eq. (11131) . any term that is symmetric under the exchange of 
indices {/t, p, p, a} does not contribute because of e pi/pa . Then, using Eq. (|Alj) . it is easy to 
see that the expansion of Eq. (1B3D to 0(£ -D / 2+1 ) in the integrand yields 


16<? 0 J e 


pupa 


f 

Jo 


dU 


-D/2+1 f d D p n _ 


(2n) 


D 


2 p 2 


x tr F vp (x)p ■ V + p ■ Vp vp {x) 


T> a + 4 T> u 


Fu P {x)p -T) + p- V7 u Jx) 


Pa 


- 2V p 7 vp {x)t> a - P pu (x)P pu (x 


(B4) 


After the momentum integrations, 


d D p 


2 p 2 


(2vr 


\D 


I PpPu I (8 tt) d / 2 \ \5 t 


(B5) 


4 


Eq. (|B4l) becomes 
16 


(8vr)^/2 


So^pupa / d££ 


—D/2+1 


tr 


F up (x)V p V a + V p V a 7 up (x) - T pv {x)F p(J [x) 


(B6) 


Finally, using Eq. (13.171) . we see that this combination identically vanishes. We infer that, 
therefore, in the pure Yang-Mills theory, 


e pi/pa G a pu (t,x)G a pa (t,x) 4 ~° (1 + 0 • gl) e puprr F“ u (x)F“ a (x) + O(t), 


(B7) 


to the one-loop order. 
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It turns out that, from very general grounds, 


wpcG a llu (t,x)G a pa {t,x ) 4 ~° e lu/pff F^ / {x)F^ (T (x) + 0(t ) 


(B8) 


holds in the pure Yang-Mills theory in all orders of perturbation theory [251 ]. To see this, one 


first notes 


M 


d t [e pupa G a pv {t,x)G a pa {t,x)] = d^Wnfax), W p {t,x) = 4 e fJ , upa D x G a Xu (t,x)G a pa (t,x), 

(B9) 

where W p (t, x) is a gauge-invariant dimension 5 axial-vector operator. This shows that 

e p , vpa G a pv (t,x)G a pc {t,x) = e pvp(T F^{x)F^ (J {x) + d ll j dt'W p (t',x). (BIO) 

Jo 

We then consider the small flow time expansion of W p (t',x) in the last term. Since there is 
no gauge-invariant axial vector of dimension < 5 in the pure Yang-Mills theory, the small 
flow-time expansion of W p (t ', x ) starts from a dimension 5 operator with an 0(t '°) coefficient 
(possibly with logarithmic corrections). This implies that the last term of Eq. (IB10D is Oft ) 
and thus Eq. (1B8[) . Our explicit one-loop calculation (1B7D is consistent with this general 
property (|B8jh as it should be. 
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